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Objectives of Project 


■ Improve simulations of soil moisture/temperature, and consequently 
boundary layer states and processes, by assimilating AMSR-E soil moisture 
estimates into a coupled land surface-mesoscale model 

■ Provide a new land surface model as an option in the Land Information 


the Land Surface 


■ Distributed land surface hydrology model 

■ Heritage: 1980’s Biosphere-Atmosphere Transfer Scheme (BATS) 

■ Can run off-line or coupled with meteorological model 

■ Flexible vertical layer configuration designed to facilitate microwave data 
assimilation 

■ Contains radiative transfer model for microwave applications 
• Described in Martinez et al. (2001), Crosson et al. 
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Ensemble Kalman Filter Data Assimilation 


1 Ne assimilate ASMR-E soil moisture observations using an Ensemble 
Kalman Filter (EnKF). Kalman filtering is a data assimilation method that 
combines a forecast (background) with observations to generate an 
improved estimate of a model variable. A Kalman Filter calculates an optimal 
weighting between the background and the observation. The EnKF uses the 
spread of the ensemble to represent the forecast error covariance. The LIS 
software includes the capability to perform EnKF data assimilation. We used 
in ensemble with 16 members generated using perturbations of 3 forcing 
variables (incident longwave and shortwave radiation, and rainfall), 14 state 
variables (14 layers of soil moisture), and 1 observation variable (AMSR-E 
soil moisture). 
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SHEELS Integration in LIS 


Wfft HP 

framework for running land surface models. 

■ We have performed off-line simulations over a Great Plains domain in LIS to 
provide initial conditions to future WRF-SHEELS coupled simulations. 

■ SHEELS ‘spin up’ has been performed off-line, forced with North American 
Land Data Assimilation System (NLDAS) data from 1/1/2002 through 6/9/2003 


Features of LIS 

■ Highly customizable at run-time, facilitating modeling experiments & 
intercomparisons 

■ Modular structure allows user to specify: 

> Land Surface Model 

► Base forcing (meteorological fields) 

► Supplemental forcing (e.g. precipitation) 

> Parameters including land cover, soil type, greenness fraction, 
topography 

' Domains of input variables may be independent. 

' Allows several tiles per grid cell to represent subgrid variability of soil type. 
1 Can run coupled with the WRF meteorological model. 


■Conically scanning passive microwave radiometer 
■Measures brightness temperatures at 6 frequencies from 6.9 to 89.0 GHz 
■Horizontally and vertically polarized radiation are measured separately at 
each frequency 

■Altitude of 705 km yields a swath 1445 km wide 

■AMSR-E/Aqua global surface soil moisture and vegetation water content are 
generated from level 2A AMSR-E brightness temperatures spatially resampled 
a nominal 25-km equal area earth grid. 

■Algorithm minimizes differences between the observed brightness 
temperatures and those generated using a forward radiative transfer model. 
■Due to extensive radio frequency interference in the 6.9 GHZ channel, 10.7 
and 18.7 GHz observations are used for soil moisture estimation. 


■ The USDA Agricultural Research Service’s (ARS) Grazinglands 

Research Laboratory has measured meteorological and hydrolsgip*. 
conditions . 

in the Little Washita Experimental Watershed - . 
network (‘ARS Micronet’) in southwestern 1} •• 

Oklahoma since 1961. •' 

■ In 1994, the ARS began monitoring the „ ’ p~ 

meteorological conditions in this watersheds. ith ' 

a network called the ARS Micronet. The Micronet •m 

now consists of 20 core stations. 

• Soil temperature and moisture are measured hourly 

at Micronet sites at 3 depths, including the 5 cm observations used 

here to validate model results. 
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SHEELS land surface model. 

Circles mark points used in 
time-depth cross sections, a 
rectangle marks the location of 
the Little Washita Micronet, 
and stars show points used 
time series at right. 
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NASA Aqua satellite with 
AMSR-E instrument 
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Isite mean. I 

I- Observations generally over-estimate soil moisture compared to the I 

Ibenchmark LIS/Staqe IV simulation and the Micronet in situ observations. B 
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lynamic range of AMSR-E observer 
moisture is small relative to that of the model, 
t correction (right) is applied to convert the 
ibservation into a model-equivalent value, 
t Cumulative Distribution Function (CDF)-matching 
technique is used here. This is similar in purpose to 
:he bias corrections usually applied to satellite 
sbservations in NWP models. 

Simulations made without the proper correction 
showed a pronounced dry bias. 
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Time series of 5 cm soil water content and soil temperature (departure from 
LIS/Stage IV simulation) at two points, and for the average of all points south of 
37N (“Mean”). 

In some cases data assimilation improves the model variables relative to the 
LIS/Stage IV benchmark. 

The two data assimilation runs generally converge. 

Effects of assimilation on soil temperatures is consistent with soil moisture 
changes - lower soil moisture results in greater diurnal range of soil 

jerature. 
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: Benchmark (LIS/Stage IV) minus DA simulation of 12-hour rainfall accumulation, 
For three model time steps just after an AMSR-E overpass. Since the DA simulation has 
in forcing set to 1.5x the Stage IV, this number is always <=0. 

Middle row: Soil water innovation (AMSR-E observed minus model soil water) for the DA 
i. There is overall good correspondence between these values and top r« 

Bottom row: Soil water 1-hr increment. This value includes the effect of data assimilation 
is well as precipitation and other physical processes during this time step. Dark red 
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:urves for different vegetation types and times of day. 

Perform more extensive validation of soil moisture and temperature simulations. 

Run WRF-SHEELS coupled model in 12-hour forecast cycles for January and July, 2003. 
Evaluate value of AMSR-E DA in estimating boundary layer states (temperature, humidity, 
vind) and surface fluxes. 

Determine landscape and hydrometeorological conditions under which assimilation is 
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